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Particle acceleration in relativistic shocks is studied analytically in the test-particle, small-angle 
scattering limit, for an arbitrary velocity-angle diffusion function D. Accurate analytic expressions 
for the spectral index s are derived using few (2 — 6) low-order moments of the shock-frame angular 
distribution. For isotropic diffusion, previous results are reproduced and justified. For anisotropic 
diffusion, s is shown to be sensitive to D, particularly downstream and at certain angles, and a 
wide range of s values is attainable. The analysis, confirmed numerically, can be used to test 
collisionless shock models and to observationally constrain D. For example, strongly forward- or 
backward-enhanced diffusion downstream is ruled out by GRB afterglow observations. 



Diffusive (Fermi) shock acceleration (DSA) is believed 
to be the mechanism responsible for the production of 
non-thermal, high-energy distributions of charged par- 
ticles in collisionless shocks in numerous, diverse astro- 
nomical systems p| . Particle acceleration is identified in 
both non-relativistic and relativistic shocks, examples of 
the latter including shocks in gamma-ray bursts (GRB's, 
where shock Lorentz factors 7 > 100 are inferred) [2|, jets 
of radio galaxies Q, active galactic nuclei Q and X-ray 
binaries (micro-quasars) 

Collisionless shocks in general, and the particle accel- 
eration involved in particular, are mediated by electro- 
magnetic (EM) waves, and are still not understood from 
first principles. No present analysis self-consistently cal- 
culates the generation of EM waves and the wave-plasma 
interactions. Instead, the particle distribution (PD) / is 
usually evolved by adopting some Ansatz for the scat- 
tering mechanism (e.g. diffusion in velocity angle) and 
neglecting wave generation and shock modification by the 
accelerated particles (the "test-particle" approximation). 
This phenomenological approach proved successful in ac- 
counting for non-relativistic shock observations. For such 
shocks, DSA predicts a (momentum p) power-law spec- 
trum, d 3 f /dp 3 oc p~ s , where s is a function of the shock 
compression ratio r, s = 3r/(r — 1) [6|. For strong shocks 
in an ideal gas of adiabatic index T — 5/3, s = r = 4 (i.e. 
p 2 d 3 f/dp 3 ocp~ 2 ), in agreement with observations. 

Analysis of GRB afterglow observations suggested that 
the ultra-relativistic shocks involved produce high-energy 
PD's with s = 4.2 ± 0.2 7]. This triggered a nu- 
merical study of test-particle DSA in such shocks, s 
calculated for a wide range of 7, various equations of 
state, and several scattering mechanisms 0, 0, 0, and 
references therein]. For isotropic, small-angle scatter- 
ing in the ultra-relativistic shock limit, where upstream 
(downstream) fluid velocities normalized to the speed of 
light c approach /3 U = 1 {(3d = 1/3}, spectral indices 
s 0; ur = 4.22 ± 0.02 were found HHHg, in accord with 
GRB observations. 
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DSA analysis is more complicated when the shock 
is relativistic, mainly because / becomes anisotropic. 
Monte-Carlo simulations H, 0, 0] and other numeri- 
cal techniques 0, 0, 0] have thus been applied to the 
problem. For isotropic, small-angle scattering, an ap- 
proximate expression for the spectrum was found |14j . 

so = (3/3„ - 2(i u pl + - (5 d ) , (1) 

and shown to agree with numerical results; in particular 
it yields so- ur = 38/9 ~ 4.22. Although s was calculated 
numerically for various scattering mechanisms [Fo| , little 
qualitative understanding of the dependence of s upon 
the scattering Ansatz has emerged. This motivates an 
analytical study of DSA in relativistic shocks, which may 
facilitate/test future self-consistent shock models. 

This letter presents the first rigorous, fully analytical 
study of DSA in relativistic shocks. The (common) as- 
sumptions made are (i) the test-particle approximation; 
(ii) small-angle scattering, given by some velocity-angle 
diffusion function D; and (iii) D is separable into an an- 
gular part times a space/momentum-dependentpart. For 
discussion and physical examples of D, see 0, H, IT^ . 
In the analysis, N angular moments of / are used (and 
higher moments neglected) to accurately calculate s even 
for small N. The moments are generalized in order to ac- 
celerate the convergence with N and to obtain s and / 
for arbitrary D. For isotropic diffusion, the analysis re- 
produces and justifies previous results such as Eq. 
good results obtained even with N — 2. For anisotropic 
diffusion, the dependence of s upon D is analyzed, and 
is numerically confirmed and complemented. It is shown 
that D can be constrained using observations, e.g. in 
GRB afterglows. The analysis works equally well for ar- 
bitrarily large 7, for which numerical methods become 
exceedingly difficult. Here we outline the main steps of 
the derivation and present its primary implications; the 
full analysis is deferred to a detailed, future paper. 

Formalism. — Consider an infinite shock front at z = 0, 
with flow in the positive z direction. Relativistic particles 
with momentum p much higher than any characteristic 
momentum in the problem are assumed to diffuse in the 
direction of momentum, jl = cos(p • z/p), according to 
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some diffusion function (DF) -D MA1 . The steady-state PD 
then satisfies the stationary transport equation |12| | 



in , - s d fifaiPu z ) d 



, (2) 



where 7 = (1 — /3 2 ) -1 / 2 is the fluid Lorentz factor and 
i G {u,d} are upstream/downstream indices, written 
henceforth only when necessary. Parameters are mea- 
sured in the fluid frame (when marked by a tilde) or 
in the shock frame (otherwise), so the Lorentz invariant 
PD f(p.,p, z) is the density in a mixed-frame phase-space. 
The absence of a characteristic momentum scale implies 
that a power-law spectrum is formed, as verified numeri- 
cally . With our assumption that is separable 
in the form Di(p,)D2(p, z), we may thus write Eq. © as 



ifi + £)8r?(£.f) = a A [(l - ?)D(ji)d& 



(3) 



where f(fi,p,z) = q(fi,f)p s , and the spatial dimension 
was rescaled as f = 7" - 1 /* D 2 (p, z)dz. The reduced DF 
D(£t) = Di(fi)/(l — jl 2 ) is constant for isotropic diffusion. 

Continuity across the shock front requires that 
fu{^u,Pu,z = 0) = fd{fid,Pd,z = 0), where upstream 
and downstream quantities are related by a Lorentz boost 
of velocity c(/3 u — /3d) /(l — PuPd)- Absence of accelerated 
particles far upstream and their diffusion far downstream 
imply that f u {z — > -00) = and fd(fid,Pd,z — > +00) = 
/oo-p^ S ' wnere foe > is constant (isotropic PD). Eq. © 
has been solved numerically under the above boundary 
conditions for specific choices of D(n) 

Angular Moments. — In the shock frame, only a small 
fraction of the particles travel nearly parallel to the flow 
(/i ~ ±1), unlike the upstream frame where a substantial 
part of / is beamed into a narrow, £l u < — 1 + 7,7 1 cone. 
Angular moments of the shock-frame PD, J*, fi n f d/i, 
thus diminish as n > increases, suggesting that the 
problem can be approximately formulated using a finite 
number of low-n moments. Next, we derive and solve 
equations for the spatial behavior of these moments; the 
boundary conditions are then used to determine s. 

We work in the shock frame and use only shock frame 
variables henceforth, incorporating the Lorentz boosts 
into the transport equations. Writing Eq. J2J) in the shock 
frame, using p = Jipi(l + Pifk) and \i = (//, + 0i)/(l + 
Pijii), we find that on each side of the shock 
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where we defined r)p s = q(p,, f)p s , D(fi) = D(ji), 
and r = j 4 f. An immediate consequence is that 



(1 — j3pL) s fi q(n, T)d/i = g = const 



(5) 



where the boundary conditions yield g u = and gd oc . 
Eq. |[5J| is an expression of particle number and energy 



conservation in the fluid frame |l4j . A more general corol- 
lary of Eq. (0} is that for any function h(fi) for which h/fi 
is finite and continuously differentiable, 



Or 



h(v)q(fJ;T)dn 
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x dp {(1 - fj 2 ) D(ji)d^ [h(n)^ 1 (1 - Pfi)- {s - 3) ] } dfi. 

The spatial dependence of an angular moment F n (t) = 
f_. fi n q(iJ,, T) d/i is given by Eq. (|BJ| with h — /j, n , if 
n > 1. Expanding the RHS integrand as a power se- 
ries around /1 = 0, we find that d T F n is given by a lin- 
ear combination of the moments F a , F a+ \ . . . Ft- Here, 
a — Max {0, n — 3}, b = n + 2 + nu, and no is defined 
as the order of the power series of D around fj, = 0, 
D(fi) oc 1 + dxfj, + g?2M 2 + ' ' ' 7 such that n^, — for 
isotropic diffusion. In matrix form, we may write 



<9 t F(t) = A 



Fo(r) 
F(r) 



(7) 



where F = (Fi, F 2 , . . .) . The (infinite) matrix A de- 
pends on 0, s and D(fi), so A„ ^ A^. The boundary 
conditions imply that {F n } are continuous across the 
shock front r = 0, vanish as t -* —00, and are finite 
for t — > +00. We show that a good, converging (in N) 
approximation to the solution is obtained using a small 
number N + 1 of moments, Fq...Fn, and neglecting 
higher order terms. This reduces F to an ./V-component 
vector and A to an N x (JV + 1) matrix. 

The spatial dependence of the zero'th moment Fq can- 
not be related to n > angular moments through Eq. 10. 
This difficulty can be circumvented by adding an addi- 
tional constraint to the system of equations. Expanding 
the integrand in Eq. J5J yields G(/3, s) ■ F(t) = g, or 



F\ - {s - 3)/3F 2 + l( s - 3)(s - 4)/3 2 F 3 + 



(8) 



independent of D(fjL). Eq. JHJ and its spatial derivative 
can be used to eliminate Fq and Fn from Eq. Q , giving 
an approximate, closed set of N = N — 1 (or less, if A is 
degenerate) ordinary differential equations (ODE's), 



9 T F(r) = A • F + gG , 



(9) 



where F = (iq, F 2 , . . . , F^) F . Here, A e and G are 
functions of A and G obtained by combining Eqs. (J7J 
and JSJ. For the relevant range of s, A is diagonalizable 
over 1Z. Thus, there exists a real matrix P G Mfj, such 
that (P _1 AP) m „ = X n S mn and the eigenvalues {A n } of 
A are real. The solution of Eq. © is then 

N 

F »( r ) ^ [^™c m e A - T -.gi 3 nm (p- 1 G) m A m 1 ] . (10) 

m— 1 

The integration constants c mjU and c mi£ j along with gd 
and s constitute 2N + 1 free parameters, as F is unique 
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only up to an overall constant. The boundary conditions 
typically impose 2N + 2 constraints, so the problem is 
over-constrained. It is natural to relax the requirement 
that -F/v be continuous across the shock, as it is relatively 
small and 8 t Fn is relatively sensitive to neglected mo- 
ments n> N. The resulting discontinuity of Fn at t = 
then provides an estimate of the approximation accuracy. 

As an example, consider the simplest non-trivial ap- 
proximation, where moments F n with n > TV = 2 are 
neglected. This is motivated, for example, by the sup- 
pression of G n [see Eq. ©], A\ tn and A 2: n for n > 2, 
if s ~ 4. Here N = 1, and the (single) ODE yields 
Fi(t) — ce Xr + gGX^ 1 . For isotropic diffusion, A = 
6/3 — (s — 2)(s — 3)(s — 4)/3 3 . In most cases of inter- 
est (where s < 5) A > 0, corresponding to exponential 
decay upstream and divergence downstream. Therefor, 
Cd must vanish, so fd is uniform in this approximation. 
This gives ct cubic ecjUcttion for s, independent of g^-, c ui 



(s - 3)/3* = |± 



(3 d - 



2(« - 



T = 



2 + (a-2)(a 



(ii) 



The real root of this equation, shown in Figure^ already 
agrees with numerical results and with Eq. Q at a ^ 5% 
level m particular it yields so:«r = 4.27. However, 
in this simple approximation Fq is not continuous across 
the shock front (it has a < 20% jump), s does not depend 
on D u (/z), and its dependence on Dd(fi) is inaccurate. 




FIG. 1: DSA for isotropic diffusion according to analytic 
approximations involving Fq . . . F2 [Eq. 1111 . heavy curves] 
or Lo-.-Ls (light curves), numerical calculations (Ref. Q, 
symbols) and Eq. Q (dotted curves). Three types of strong 
shocks (see ^]) are examined: using the Jiittner-Synge equa- 
tion of state (solid curves/crosses), assuming a fixed adiabatic 
index r = 4/3 (dashed curves/x-marks), and assuming a rel- 
ativistic gas where /3 u (3d = 1/3 (dash-dotted curves/circles). 

The convergence of the approximation with N depends 
on D, more moments required in general as njj and |d„| 
increase. Due to degeneracies in A, N = 2, 3, 4, 5, 6, . . . 



correspond to N = 1,1,2,3,3,..., and N > 6 is re- 
quired in order to improve the N = 2 approximation, e.g. 
so-,ur(N = 6) = 4.24. The convergence is significantly ac- 
celerated if the moments are properly generalized. 

Generalized Moments. — The analysis can be general- 
ized in many ways by exploiting the freedom in the defi- 
nition of the angular moments. For example, if we define 
Fn{r) = j\ q{ii,T)w n {n)J n {n) dfj,, where {/„} is an or- 
thonormal basis in the interval — 1 < fj, < 1 with weight 
functions {w n }, then the moments {F n } are simply the 
expansion coefficients of a in terms of the basis {/«.}, i.e. 
q(fx,r) = S n F„(r)/ n (/i) With this choice of 

the analysis provides direct information about the PD /. 
Additional constraints on / [e.g. [l4| may thus be used 
to estimate/improve the accuracy of the approximation. 
The analysis is mostly unchanged; it remains analytically 
tractable as long as A, G can be determined analytically, 
although the analogue of Eq. (|llfl is in general no longer 
a polynomial, becoming a transcendental equation in s. 

As an example, let /„(/i) = \{2n + l)/2] 1 ^ 2 P n (^) and 
w n = 1 , where P n (fj,) is the Legendre polynomial of or- 
der n. The moments {F n } thus equal the coefficients 
{L n } of the Legendre series of q. Here we cannot take 
N = 2, because the corresponding eigenvalue A is neg- 
ative. For N = 3, 4, 5, . . . we find N = 2,3,4,..., with 
N u = 1,2,2,... (Nd — 1,1,2,...) non-vanishing expo- 
nents upstream (downstream). The resulting spectrum 
for isotropic diffusion is accurate to ~ 5%, 2%, 1%, . . .; 
in particular So ; « r = 4.16,4.21,4.22,.... These results 
illustrate the rapid convergence of the Legendre-based 
analysis. The method is equally applicable for arbitrarily 
high 7, whereas numerical methods become exceedingly 
difficult as the ultra-relativistic limit is approached. 

Isotropic (anisotropic) diffusion results for s with N = 
5 Legendre moments are shown in Figure ^ (Figure |5J). 
For isotropic diffusion, previous results are reproduced, 
lending rigor support for Eq. (Q. The lowest order mo- 
ments converge quickly, giving a rough estimate of g(/x, r) 
as a sum of components exponentially decaying in |r| and, 
in the downstream, also a uniform part. The analytically 
calculated s, q compare favorably with numerical results, 
which we calculate in the eigenfunction method of 0. 

Anisotropic diffusion. — Qualitatively good results for 
d n <3 7^ 0; d n >4 = are already obtained with N = 3 
Legendre-based moments. More complicated forms of the 
DF, involving additional d n ^ terms, require progres- 
sively larger N. Figure |21 shows that large As = s — s 
values can be obtained for anisotropic diffusion in rela- 
tivistic shocks. The spectrum is more sensitive to Dd 
than it is to D u , by a factor of a few. It is more sensitive 
to di than to d2] \As(D = l + d n n n )\ is roughly monoton- 
ically decreasing with n. A linear, forward- (backward-) 
enhanced downstream DF, Dd = 1 + di/i, yields s > 4.3 
(s < 4.0) in a 7 > 10 shock, if d x > 0.6 (d x < -0.85). 
Some forward- (backward-) enhanced choices of D lead to 
more extreme spectra, e.g. s ur > 4.6 (s ur < 3.7). Con- 
straints on Dd may thus be imposed using observations 
that rule out such s values, e.g. in GRB afterglows. 
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FIG. 2: DSA for anisotropic diffusion. The spectrum is cal- 
culated numerically /with L0 . . . L5 (symbols/curves) for D = 
f + d n /j, n in a 7 = 10, strong shock with the Jiittner-Synge 
equation of state. Plotted is As = s — so vs. d n for n = 1 (cir- 
cles/solid curves) and n = 2 (squares/dashed curves), down- 
stream/upstream (filled/open symbols, heavy/light curves). 
Inset: local deviations from an isotropic DF down/upstream 
(heavy /light curves): D — 1 + do ©(too — \fi— /j,q\) (solid curves) 
and D = 1 + do exp[— (fj, — /io) 2 /2wq] (dashed curves), where 
(do, wo) = (0.1,0.05) and O is the Heaviside step function. 



In order to elucidate the role of D, we numerically cal- 
culate s for small deviations AD, localized around some 
angle /iq, from an otherwise isotropic DF. The parame- 
ter As/AAd, where AAd = J jADdfj,, depends only 
weakly on the exact form of such AD. Figure El (inset) 
shows that s is sensitive to fj,o, especially downstream. 
In general, the spectrum hardens when Dd is enhanced 
at [i < fi c ~ 0.4 or diminished at /i > /i c , and vice 
versa. Angles for which /1 < /i c (/1 > /i c ) roughly cor- 
respond to negative (positive) flux in momentum space, 
cx — D(/i)(l — fi 2 )d fl q, qualitatively accounting for this 
behavior. One way to see this is to note that the spec- 
trum hardens when the average fractional energy gain 
/e (of a particle crossing, say, to the downstream and 
returning upstream) or the return probability (to the 
upstream) P ret increase; s = 3 + ln(l/P ret )/ln/£; [To| . 
For example, enhanced diffusion in angles where j M < 
shifts qd towards the upstream, raising both P ret and Je- 
The above trend is roughly reversed upstream, but As is 
smaller there and more sensitive to details of AD U . 
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